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Abstract 

A  method  of  determining  pressure  changes  in  a  two-phase  water  system  is 
presented.  Two  applications  are  given  which  pertain  to  the  design  of  nuclear 
reactor  system.s.  The  method  utilizes  informiation  obtained  from  simple  calcu¬ 
lations  to  enter  a  chart  from  which  the  pressure  can  be  readily  determined. 

In  addition,  an  anal3rtical  approach  is  taken  and  an  equation  derived  which 
describes  pressure  changes  also,  but  less  accurately  than  the  graphical 
method . 

Nomenclature 


a  -  constant,  defined  as  used 

b  -  constant,  defined  as  used 

c  -  constant,  defined  as  used 

Cy  -  specific  heat  of  air  at  constant  volume,  Btu/lb-°F 

d  -  constant,  defined  as  used 

h  -  enthalpy,  Btu/lb 

lb  -  pound  mass 

m  -  constant,  defined  as  used 

n  -  constant,  defined  as  used 

P  -  pressure,  absolute 

Q  -  energy  added  by  heat  transfer,  Btu 

t  -  temperature, 

u  -.internal  energy,  Btu/lb 

V  -specific  volume,  ftVlb 

V  -  internal  volume  of  a  vessel,  ft-^ 

w  -  mass  ,  lb 

^  -  pound  force 

Subscripts 

amb  -  initial  ambient  value 

av  -  average 

f  -  refers  to  liquid  phase 

fg  -  property  of  vaporization 

g  -  refers  to  vapor  phase 

in  -  m,ass  or  energy  entering 

out  -  mass  or  energy  leaving 

P  -  at  constant  pressure 

-  at  constant  intern^,,  energy 

-  at  constai  t  specif?!^  volume 

4  1  -  initial  state 

,  2  -  final  state 
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Graphical  Method 


It  is  well  known  that  if  two  intensive  thermodynamic  properties  of  a  p:re 
substance  in  equilibrium  are  known,  then  all  the  other  like  properties  are  fixed 
and  can  be  found  as  functions  of  the  known  properties.  The  pure  substance  ex¬ 
amined  herein  is  water  in  two  phases  ( 15 quid  and  vapor)  which  are  in  equilibria. 
The  phases  are  contained  in  a  rigid ,  non-adiabatic  vessel  which  is  connected  to 
sources  and  sinks  of  the  same  substance  not  necessarily  saturated.  The  initial 
amounts  of  liquid  and  vapor  are  known  as  well  as  the  initial  pressure.  The  as¬ 
sumption  is  made  that  electricity,  magnetism,  capillarity,  changes  in  elevation, 
and  kinetic  energy  are  negligible.  Consider  now  that  some  fluid  leaves  the 
vessel  (composed  of  either  or  both  phases)  and  is  replaced  in  part  or  in  whole 
by  a  mass  of  liquid  or  vapor  in  an  Indefinite  amount  of  time.  Referring  to 
Figure  1,  and  treating  the  vessel  as  a  control  volume,  the  following  balances 
may  be  written! 


Energy  Balance 

'^f2"f2  V^g2  =  '^fl^^fl  ^  W  ^  "  '^out^out  +  ^in 

Mass  Balance 

'"f2  +  V  =  '"fl  '■'gl  +  ’^in  -  ^out  ^2) 

VqXymq  ”Ba3,anqe" 

V'"g2  ' 

One  method  of  solving  for  the  final  pressure  is  by  the  cut-and-try  process. 
A  final  pressure  is  assumed  and  the  values  of  internal  energy  and  specific  vol¬ 
ume  are  gotten  from  a  table  of  properties  such  as  reference  (l).  Equations  (l) 
■and  (3)  can  then  be  solved  for  wp2  ^-f^d  ^^2'  Souation  (2)  is  then  used  as  a 
check  on  the  correctness  of  the  assumed  pressure. 


This  process  can  be  reduced  to  a  straight-forward  solution  by  utilizing 
the  definition  of  steam  quality,  "x”. 


X  =  Wg/(w^  +  Wg)  =  w^/w 

From  this  definition,  the  following  is  derived: 

(4) 

Wg=XW 

(4a) 

wp  =  (1  -  x)w 

(4b) 

Furthermore , 

U  =  (1  -  x)u^  +  xu^ 

(5) 

and, 

V  =  (1  -  x)vp  •+  XVg 

(6) 
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utilizing  equations  (4n),  (4b),  (5)  and  (6),  equation  (l)  can  be  rewritten  as: 

^2^2  ~  '“^in'^in  “  '^out^ovt  '^in 

Similarly,  equation  (2)  can  be  written: 


W2  =  Wi  +  W.^  -  w^ut 

Equations  (V)  and  (S)  can  be  solved  and  combined  to  find  ^2' 
^2  -  W2^i2/'^2 


(?) 

(9) 


Also: 


V2  =  V/w2 

Vie  now  have  two  intensive  properties  and  know  that: 

P2  =  P(u2,V2) 

By  entering  the  u-v  chart  (Figure  2)  with  the  values  of  U2  and  V2,  the  pressure 
can  be  determined. 


Example  I  -  Pressurizer  Outsurtre . 

To  illustrate  the  use  of  the  foregoing  method,  let  us  assume  that  we  wish 
to  know  how  fn.r  the  pressure  in  a  pressurized  water  reactor  system  will  fall 
during  an  increasing  power  demand  transient.  The  mass  of  liquid  which  leaves 
the  pressurizer  will  be  equal  to  the  mass  required  to  make  up  .the  contraction 
in  the  main  coolant  system.  During  rapid  transients,  the  effect  of  pressurizer 
heaters  is  negligible  and  equation  (V)  becomes: 

W2U2  =  w^u^  -  (12) 

where  hf  ,0^  is  "the  average  enthalpy  of  the  pressurizer  liquid  during  outsurge. 
This  average  value  can  usually  be  approximated  quite  closely,  bi’t  if  very  accu¬ 
rate  results  are  desired,  two  or  more  iterations  will  be  necessary.  Equation 
(?)  becomes: 

'•^2  =  ^1  -  Vrt 

Dividing  (l2)  by  (l3)  gives  U2.  Dividing  the  volum.e  of  the  pressurizer  by  the 
result  of  (13)  will  yield  V2-  The  final  pressure  is  found  by  entering  Figure  2 
with  these  quantities. 

Example  II  -  Vauor  Container 

The  pressurized  water  reactor  system  of  the  present  day  is  always  designed 
to  be  installed  in  an  enclosure  which  will  contain  the  vapor  resulting  from  a 
loss  of  coolant  accident.  The  usual  design  considers  that  all  the  fluid  in  the 
primary  system  is  lost  as  well  as  the  fluid  in  one  or  more  of  the  associated 
steam  generators.  In  general,  there  are  five  fluids  to  be  considered,  namely: 

1.  Primary  coolant  in  the  main  coolant  system  and  auxiliary  prim.ary 
systems . 

2.  Primary  liquid  in  the  pressurizer. 
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3.  Primary  vapor  in  the  pressurizer. 

4.  Secondary  liquid  in  the  steam  generator(s) . 

5.  Secondary  vapor  in  the  steam,  generator(s)  . 

It  is  usually  assumed  that: 

1.  No  condensation  of  vapor  takes  place  on  the  walls  of  the  containm.ent 
vessel. 

2.  No  energy  is  added  to  the  vapor  from  hot  m.etal  in  the  plant. 

3.  The  resulting  liquid  and  vapor  phases  are  in  equilibrium  with  each 
other  and  with  the  air  in  the  vessel. 

4-  All  the  fluid  considered  is  released  instantaneously. 

5.  No  vapor  or  air  escapes  from,  the  vessel. 

The  energy  balance  for  the  resulting  liquid-vapor-air  system,  may  be  written 
as  follows: 


”2^2  "  fluid)  -  w^^^c.^^(t2  -  (l4) 

where  t2  is  the  saturation  temperature  corresponding  to  the  final  partial  pressure 
due  to  steam.  The  mass  balance  is: 


W2  =  Su  (of  each  fluid) 


(15) 


It  is  necessary  to  start  with  a  vapor  container  of  given  design  and  compute 
its  void  volume,  V^.  The  final  internal  energy  and  specific  volume  are  compacted 
as  before  and  the  final  pressure  of  the  vapor  obtained  from  Figure  2.  To  this 
must  be  added  the  partial  pressure  due  to  the  heated  air; 

^2, air  “  14*'7(t2  +  460)/(tg^nib  4^0) 

The  total  absolute  pressure  is  then  the  sum  of  the  partial  pressures  due  to  vapor 
and  air. 

It  would  seem,  from  equation  (I4)  that  the  final  pressure  must  be  known  before¬ 
hand  in  order  to  have  a  value  of  t2  and  that  therefore  a  cut-and-try  process  m.ust 
be  used.  This  is  true  only  if  precise  answers  are  desired.  Practically,  vapor 
containers  are  designed  for  low  pressures  of  the  order  of  50  to  100  psi  so  that 
a  good  approximation  may  be  made  to  t2.  Furthermore,  the  second  term  on  the  right- 
hand  side  of  (14)  is  very  small  compared  to  the  energy  in  the  plant  fluids.  There¬ 
fore  a  rather  large  error  in  estimating  the  final  temperature  would  have  to  be  made 
in  order  to  significantly  affect  the  result. 

Analytical  Approach 

Because  pressure  is  a  point  function,  its  differential  is  exact.  Ne  can 
therefore  write  from  (ll)  the  following: 

dP  =  (^P/^u).y.du  +  (^P/^v).ydv  (!.'■;) 
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From  the  calculus  we  know  that; 


(^Pu)^(V^)p(<=^v/»^  =  -1 


0?/)n)^  =  -(Vc)tOp(‘^PA)t, 

Also,  for  saturated  states: 

QvA»)p  =  Vfg/u^g 

Utilizing  equations  (lU)  and  (l9) ,  equation  (l6)  may  be  written: 

dP  =  (dv  -  j^v^g/ufg  j  du)  ^20) 

Integrating  (20)  from  state  1  to  state  2  and  assuming  that  0  pA)^  is  a  constant! 

aP  =  (av  -  [vf.g/ufg]Au)PpA)y  ^21) 

Practically  speaking,  equation  (2l)  is  of  limited  usefulness  because  of  the 
assumption  that  (<^P/cV)^  is  constant.  An  examination  of  Fi^ire  3  will 
that  this  property  varies  markedly  even  for  relatively  small  changes  in  pressure. 
It  is  possible  however,  to  linearize  over  useful  ranges  by  writing; 

=  a  +  bP  (22) 

Substituting  (22)  into  (20)  and  integrating  from  states  1  to  2: 


=  (a/b  +  P^)  -  a/b 


The  properties  Vf  and  Uf„  are  to  be  evaluated  at  the  initial  pressure. 
(^/c)v)y  is  to  be  evalxiated  at  the  final  Internal  energy.  For  an  explanation 
of  this,  the  reader  is  referred  to  Figure  4  and  the  following  paragraph. 

„  Examination  of  equation  (2l)  shows  that  the  term  -(v^g/u^  )(au)  competes 
v'  -  v,  and,  since  pressure  is  constant,  defines  state  point  1  .  M\iltiplying 
(vj  -  v^^)  by  (5pA)y  results  in  (Pj^  -  P')  ,  and  defines  state  point  1*  since 
internal  energy  is  held  constant  at  U2.  The  term  Av(oP/:)v)■^.  gives  P  P2- 
The  sum  of  these  two  differences  is  then  Pj  -  P2' 

Observing  that  (’^P/cV).,,  vs  P  is  nearly  linear  on  the  log-log  plot  of 
Figure  3  suggests  that  a  more  general  equation  could  be  vrritten  by  letting 
this  property  assume  the  form: 

=  ap”  (24) 

It  is  evident  from  Figure  3  that  "a"  and  ”b"  are  themselves  functions  of  the 
internal  energy.  These  can  be  represented,  to  a  good  approximation,  by: 

_  ,  (25) 

a  =  mu  +  c 
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b  =  nu  +  d 


(26) 


and , 


Substituting  (24),  (25)  and  (26)  into  ( 20),  evaluating  m,  n,  c,  and  d  from  Figure 
3,  and  integrating.,  the  following  equation  is  obtained: 


^2  = 


-1.4  +  4.5x10"'^2 


+  (-1.4  +  U.5ylO~\:^)/{e.O  -  5.:5x30"^) 


(av 


l/(-1.4  +  4.5x10 


(27) 


Equation  (27)  is  rather  cumbersome  for.nanmal  calculation,  but  could  be  use¬ 
ful  in  computer  solutions  where  charts  obviously  cannot  be  used.  The  constants 
used  in  (27)  were  computed  on  the  basis  of  a  best  fit  for  Fimre  3  over  its  entire 
range.  Better  approximations  could  be  made  by  limiting  the  range  to  one  of  inter¬ 
est  for  a  particular  problem.  Figures  2  and  3  were  drawn  on  the  basis  of  proper¬ 
ties  given  in  reference  1. 
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